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ABSTRACT 

Investigating the spin parameter distribution of subhaloes in two high resolution isolated halo simulations, re¬ 
cent work by Onions et al. suggested that typical subhalo spins are consistently lower than the spin distribution 
found for field haloes. To further examine this puzzle, we have analyzed simulations of a cosmological volume 
with sufficient resolution to resolve a significant subhalo population. We confirm the result of Onions et al. and 
show that the typical spin of a subhalo decreases with decreasing mass and increasing proximity to the host 
halo center. We interpret this as the growing influence of tidal stripping in removing the outer layers, and hence 
the higher angular momentum particles, of the subhaloes as they move within the host potential. Investigating 
the redshift dependence of this effect, we find that the typical subhalo spin is smaller with decreasing redshift. 
This indicates a temporal evolution as expected in the tidal stripping scenario. 

Subject headings: methods: numerical - galaxies: haloes- galaxies: evolution - cosmology: theory - dark 
matter 


1. INTRODUCTION 

In the standard model of structure formation, the rotation 
velocities of disc galaxies are correlated with the spin prop 


erties of their surrounding dark matter haloes (Fall & Efs- 
tathiou|[l980 1. The simplest model explains this correlation 


via angular momentum conservation and assuming baryons 
and dark matter initially share the same specific angular mo¬ 
mentum distribution ( |Mestel|1963[ ). Even though subsequent 
models pa int a more complex pict ure, this link continues to 
exist (e.g. Dalcanton et al. 1 1997 Mo et al.|1998| Navarro & 


ISteinmetz 2000||Abadi et al.|2063||Bett et al.|2010) . Conse 


quently, the halo spin is an important parameter in many semi- 


analytic models of galaxy formation (I 

Cauffmann et al. 

1993 

1997llFrenk et al. 1997; Cole et al. 20C 

)0; Eenson et al. 

2001 

Bower et al. 2006; Croton et al. 2006 

De Lucia & Blaizot 

2007 Bertone et al. 2007; Font et al 

2008, Benson 2012) 


vidual dark matter haloes in cosmological simulations (Pee- 


bles 1969 

Bullock et al.| 2001; Hetznecker & Burkert 

2006 

Bett et al. 

2007 Maccio et al .12007 Gottlober & Yepes 

2007 

Knebe & 

Powerj20081 

Antonuccio-Delogu et al. 2010; 

Wang 

et al. 201 If iLacerna & 

c Padilla][20l2| Trowland et al. 

2013; 


Bryan et al. 


Due to a lack of resolution in previous generations of large 
cosmological simulations, subhalo spins have not been thor¬ 
oughly investigated so far, despite their application within 
cu rrent semi-analy t ic mo dels ( |Guo et al.|2011| ). Initial work 
by|Lee & Lemson|(|2013) analyzed the spins of the two most 
massive substructures of Local Group like systems in the 
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Millennium-II simulation (Boylan-Kolchin et al7j[2009|) and 
revealed possible consequences for the application of subhalo 
spins to near-field cosmology. 

Onions et al. ( |2013| ) investigated the spin distribution of 
subhaloesln two high resolution simulations of a Milky Way- 


like halo (the Aquarius simulation ( 

Springel et al. 2008) and 

the GHALO simulation (Stadel et al. 

2009 

)) analyzed by a va- 
that subhalo spins 

riety of subhalo finders. They suggested 


are significantly offset to lower values than those seen in typi- 
cal distribution f unctions fitted to haloes (Bullock et al.|200T| 
|Bett et aT7||2007| >. As this result is independent of the sub¬ 
halo finder used, it suggests that this is a true physical effect. 
This could not be investigated further because their resimula¬ 
tion did not contain a l arge fie ld halo population. Excluding 
subhaloes, |Colfn et al.| (|2004J) found the spin parameter dis¬ 
tribution of isolated dwarf dark matter haloes to be perfectly 


: perfectly 

consistent with that of larger haloes. This suggests that the 
consistently lower spin of substructure is not due to the gen¬ 
erally smaller mass of subhaloes, but is more likely related 
to tidal stripping of high angular momentum material. On 
the other hand, this offset could also be due to differences be¬ 
tween the Aquarius and GHALO simulations and those used by 
Bullock to define the field relation. To answer this question 
we require a single simulation that simultaneously includes 
both a significant subhalo and field halo population. 

In this work we use purpose built simulations, specifically 
designed to contain both a field and subhalo population, to in¬ 
vestigate the difference in spin distribution functions between 
subhaloes and haloes. In Section|2j we present these simula¬ 
tions and the corresponding (sub-)halo catalogues. The differ¬ 
ent theoretical models of dimensionless spin parameters are 
described in Section [3] while our results are summarized in 
Section|4] We discuss our work and conclude in Section^ 

2. SIMULATION DATA 


As we require our halo and subhalo masses to span a 
wide dynamic range (10 8 < M < 10 15 Mg) , we have 
run four dark matter only comoving cosmological boxes 
containing 512 3 particles, with linear sizes of 8/i -1 Mpc, 
20/i _1 Mpc, 50/i -1 Mpc, 100/i -1 Mpc respectively (here¬ 
after BoxA, BoxB, BoxC, BoxD). The softening lengths are 
chosen to be 4% of the mean separation between particles. 
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This set of simulations can both sufficiently resolve subhalo 
spins (at least 300 particles per subhalo; Bett et al. 2007) 


and have significant statistics for haloes (c.f. Table [T i. We 
also ran two simulations with the same parameters as those 
used for the BoxA simulation except for the gravitational soft¬ 
ening parameter. BoxA_Sl has a smaller softening length 
whereas BoxA_S2 has a larger softening length. We also ran 
a low resolution simulation containing 256 3 particles and the 
same linear size, 8/i -1 Mpc as BoxA which we designate 
BoxLo. The mass resolution of BoxA (2.6 x 10 5 /i _1 Mq 
per particle) is very close to that of the Aquarius-A simu¬ 
lation at level 4 (Springel et al. 2QQ8j ) which had a parti¬ 
cle mass of 2.7 x 10° h 1 Mq in the High resolution region. 
This is roughly three times better mass resolution than the 
Millenni um-II sim ulation ( m p = 6.9 x 10 6 used 

by Lee & Lemson ( |2013[ ). The cosmology was chosen to 
be the same as in the Aquarius simulation, i.e. ACDM with 
Qm = 0.25, n A = 0.75, <7g = 0.9, n s = 1, h = 0.73. 
Initial conditions were generated at z = 127 by the code N- 
GenlC using the Zel’dovich approximation (written by Volker 
Springel) to linearly evolve positions from an initially glass¬ 
like state. This was then evolved to the present day using 
gadget-2 (Springel|2005| ). 


All simulations except BoxLo were analyzed with the (sub- 


)halo finding code SUBFIND (Springel et al. 

2001), AHF (Gill 

et al. 

2004 

Knollmann & Knebe 2009) and 

HBT (Han et al. 


2012). BoxLo was only analyzed with SUB FIND . A summary 
of our simulations is given in Table |T| and details of all the 
halo finding algorithms we have used and a di scus sion of their 
relative merits can be found in [Knebe et al.[ ( {20T3] ). 

3. THEORY 

The dimensionless spin parameter indicates how much a 
collection of particles is supported by the angular momentum 
against gravitational collapse assuming gravitational equilib¬ 
rium, where a negligible spin parameter represents minimal 
support, while the value for a completely supported system 
depends on the chosen parametrisati on. There are t wo stan¬ 
dard parame trisations defined by Peebles ( |1969| ) and Bullock 
et aL|(200ljl, respectively, that we describe in the following 


two sections. 

|Hetznecker & Burkert| ( |2006| ) showed that Bullock’s 
parametrisation is less dependent on redshift evolution than 
Peebles’ parametrisation. This is due to it being more robust 
to variations in the position of the structure’s outer radius and 
therefore not as strongly affected by the many minor mergers 
over a halo’s merging history. Therefore, the two descriptions 
are not readily interchangeable and results need to be com¬ 
pared using the same parameter. 

3.1. Peebles Spin Parameter 

Peebles (1969) proposed to parameterise the (sub-)halo spin 
in the following way: 


A = 


jyw\ 

GM 5/2 ’ 


( 1 ) 


where J is total angular momentum, E the energy and M the 
mass of the (sub-)structure. With this choice a value of A ~ 
0.4 represents a purely rotationally supported object (Frenk & 
|White|2012| ). 

Applying this parametrisation Bett et al. (2007) determined 
the spin distribution of haloes in the Millennium Simulation 
(Springel et al.|2005). The Millennium Simulation has a mass 


resolution of m p = 8.6 x 10 s Ii~ 1 Mq and therefore contains 
very few subhaloes. The vast majo rity of the objects in the 
TREEclean catalogue of Bett et al. (2007| ) are haloes rather 
than subhaloes. The probability density function of log A they 
found to fit the distribution used the following parametrisa¬ 
tion; 


P(log A) = A ( ) exp 


3/a 


where A is given by. 


A = 3 In 10 


T(a)’ 


( 2 ) 


(3) 


and T(a) is the gamma function. They found Ao = 0.04326 
and a = 2.509 best fit the distribution of halo spins. 

3.2. Bullock Spin Parameter 

Bullock et al. ( 2001f proposed a dimensionless spin param¬ 
eter of the form: 

A' = -jJ- -, (4) 

y/2 MRV 

where J is the angular momentum within a virilized sphere 
with radius R and mass M, and V is the circular velocity 
at the virial radius ( V 2 = GM/R). They also proposed a 
parametrisation of the probability density function based on 
|Barnes & Efstathiou j ( j 1987) , 


P(A') = 


1 

X'\/2ttc 


exp - 


ln 2 (A7A') 

2er 2 


(5) 


They found the best fit for haloes is given by X' 0 = 0.035 and 
er = 0.5. 

4. RESULTS 

Note that in this section we will show results for the dis¬ 
crete, normalized derivative of the spin distribution function 
AN(< logA)/A log X/Ntot and A N(< A')/AX'/N totaU 
while the fitted functions are for the continuous probability 
density function P(logA) = dN(< log A)/dlog X/N to tai 
and -P(A') = dN(< X')/dX'/N to tai, respectively. For Pee¬ 
bles spin, we set the bin width to be A log A = (log A max — 
logA min )/100. For Bullock spin, the bin width is AA' = 
(A ' max - A^J/100. 

4.1. Halo finding code, softening, and resolution test 

In this section, we first test whether simulation data set 
and the specific choice of substructure finding code will af¬ 
fect the derived spin of (sub-)haloes. Maccio et al. (2008) 
has tested halo spins with different cosmological parame¬ 
ters. They found that the spin distributions of haloes is essen¬ 
tially independent of cosmology, at least for changes between 
WMAP1.WMAP3 and WMAP5. We choose not to confirm 
this result here. BoxA, BoxAJSl and BoxA S2 are used to 
compare different force resolutions. BoxA and BoxLo are 
used to compare different mass resolutions. To ensure reliable 
properties are recovered only (sub-)haloes with more than 300 
particles are selected throughout this work. While calculating 
the spin of haloes, all their substructures are removed from 
them. Spin distributions are fitted by Equation [2] and Equa¬ 
tion^ and the fitting parameters are listed in Table[2] 
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Table 1 

Summary of simulation properties, halo and subhalo counts for the halo finder indicated. 


Name 

Box size 

Particle mass 

Force softening 

Nhalo, SUBFIND 

Nhalo, HBT 

Nhalo, AHF 

N S ub, SUBFIND 

N S ub, HBT 

Nsub,AHF 


h~ 1 Mpc 

h 1 Mg 

/i —1 kpc 

> 300 

> 300 

> 300 

> 300 

> 300 

> 300 

BoxLo 

8 

2.1 x 10 6 

1.25 

1136 

- 

- 

213 

- 

- 

BoxA.Sl 

8 

2.6 x 10 5 

0.04 

6589 

6698 

6775 

1934 

2169 

1460 

BoxA 

8 

2.6 x 10 5 

0.63 

6651 

6587 

6798 

1651 

1899 

1216 

BoxA_S2 

8 

2.6 x 10 5 

1.25 

6585 

6476 

6529 

1388 

1618 

944 

BoxB 

20 

4.1 x 10 6 

1.56 

8923 

8785 

9139 

2111 

2494 

1302 

BoxC 

50 

6.5 x 10 7 

3.91 

12791 

12533 

12874 

2687 

3325 

1597 

BoxD 

too 

5.2 x 10 s 

7.81 

17562 

17053 

16901 

3072 

3949 

1737 





> 2400 

> 2400 

> 2400 

> 2400 

> 2400 

> 2400 

BoxA* 




1132 

- 

- 

215 

- 

- 


This higher particle number threshold is used to compare with BoxLo over the same halo mass range. 


Table 2 

Parameters for the spin distribution with different substructure finding codes and force resolution. 


Peebles Spin 

haloes 

Ao 

a 

subhaloes 

Ao ol 

SUBFIND 

BoxA_Sl 

BoxA 

BoxA S2 

0.0398 ± 0.00019 
0.0371 ± 0.00021 
0.0364 ± 0.00026 

2.54 ±0.039 
2.59 ±0.047 
2.59 ±0.058 

0.0237 ± 0.00032 
0.0254 ± 0.00049 
0.0298 ± 0.00036 

3.56 ±0.13 
2.92 ±0.17 
2.60 ± 0.10 

HBT 

BoxA.Sl 

BoxA 

BoxA S2 

0.0390 ± 0.00024 
0.0366 ± 0.00019 
0.0356 ± 0.00021 

2.42 ± 0.049 
2.48 ± 0.042 
2.37 ± 0.046 

0.0260 ± 0.00029 
0.0277 ± 0.00027 
0.0322 ± 0.00033 

2.89 ±0.10 
2.68 ±0.08 
2.52 ±0.08 

AHF 

BoxA.Sl 

BoxA 

BoxA_S2 

0.0380 ± 0.00024 
0.0369 ± 0.00026 
0.0367 ± 0.00023 

2.82 ±0.053 
2.76 ±0.059 
2.61 ± 0.051 

0.0303 ± 0.00044 
0.0343 ± 0.00045 
0.0384 ± 0.00074 

2.53 ±0.12 
2.84 ±0.11 
3.09 ±0.17 

Bullock Spin 

haloes 

<T 

subhaloes 

A' a 

SUBFIND 

BoxA.Sl 

BoxA 
BoxA S2 

0.0308 ± 0.00027 
0.0310 ± 0.00025 
0.0320 ± 0.00022 

0.655 ±0.007 
0.629 ±0.007 
0.615 ±0.006 

0.0123 ±0.00017 
0.0167 ± 0.00019 
0.0224 ± 0.00028 

0.827 ± 0.009 
0.664 ± 0.009 
0.629 ±0.010 

HBT 

BoxA.Sl 

BoxA 

BoxA S2 

0.0303 ± 0.00026 
0.0309 ± 0.00026 
0.0320 ± 0.00025 

0.638 ± 0.007 
0.626 ±0.007 
0.605 ± 0.007 

0.0150 ±0.00013 
0.0198 ±0.00018 
0.0262 ± 0.00031 

0.754 ± 0.007 
0.674 ± 0.008 
0.644 ± 0.010 

AHF 

BoxA.Sl 

BoxA 

BoxA_S2 

0.0278 ± 0.00025 
0.0283 ± 0.00019 
0.0290 ± 0.00022 

0.677 ± 0.008 
0.655 ±0.006 
0.613 ±0.006 

0.0178 ± 0.00028 
0.0248 ± 0.00046 
0.0323 ± 0.00059 

0.768 ±0.013 
0.761 ± 0.016 
0.779 ±0.015 


As Table [2] shows, the recovered spin properties of haloes 
are largely independent of the choice of gravitational soften¬ 
ing. For subhaloes there is a slight trend for the Peebles spin 
parameter to increase as the softening is increased but this ef¬ 
fect is only barely resolved. Such a trend would be expected 
as a larger softening will produce a shallower core potential, 
lowering slightly the central kinetic energy and altering the 
energetics and angular momentum profile, thus affecting the 
spin parameter. For the Peebles measure (Equation IT]), the 
change in energy is outweighed by the change in angular mo¬ 
mentum but the two effects counteract each other. For the 
Bullock spin parameter (Equation |4], only angular momen¬ 
tum has an affect on the spins. So we find that the Bullock 
spin of subhaloes is more sensitive to the softening, as shown 
in Table [2] 

Table|2]also shows that, contrary to|Onions et al.|(J2013jl, the 
three halo finding methods do not recover consistent spin pa¬ 
rameters. While they all agree on the halo spins, AHF recovers 
significantly larger spins on average for the subhalo popula¬ 
tion than either SUBFIND or HBT which are consistent with 
each other. The subhalo spins for AHF are broadly consistent 


with the field population, particularly for larger values of the 
gravitational softening. This is also discrepant with Onions 


et ak (12013) who found lower spin parameter values for their 
subhaloes - This result is due to the failure of AHF to resolve a 
significant fraction of subhaloes within the simulations. The 
subhalo numbers given in Table [I] indicate that around 36% 
of the subhaloes containing 300 or more particles in BoxA 
are missed by AHF. Difficulties for AHF in resolving sub¬ 
structures where the density contrast between the subhalo and 
the main halo is expected to be small have also been reported 
elsewhere ( jAvila et al.|2014] ). Further evidence for this issue 
is the rising incidence of missing substructures as the box size 
is increased evidenced in Table [l] for the largest box (BoxD), 
AHF misses 71% of the subhaloes found by SUBFIND. AHF 
is missing small subhaloes in the outskirts of the host halo 
and (as we shall demonstrate later) these small subhaloes are 
precisely the ones with the lowest spin parameters. 

This naturally produces a population of subhaloes with 
higher average spin parameter for AHF. As HBT and SUBFIND 
produce consistent results and AHF fails to recover the com¬ 
plete subhalo population we choose to concentrate our analy- 
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sis on SUBFIND for the remainder of this paper. 

In Table [3 we test the influence of mass resolution. In 
the brackets after simulation name, we note the particle num¬ 
ber threshold chosen. This is set in order to match (sub)halo 
masses between BoxA and BoxLo, the lower resolution ver¬ 
sion of this simulation. This ensures that the halo and subhalo 
catalogues for BoxLo(lVp > 300), BoxA(7V p > 2400) and 
BoxA S2(A p > 2400) have the same mass range. The gravi¬ 
tational softening lengths for both BoxLo and BoxA were set 
to 4% of mean particle separation. BoxLo and BoxA_S2 have 
the same absolute softening length (1.25/i -1 kpc). The results 
show that mass resolution has almost no effect on spin dis¬ 
tribution. The first three rows of each sub-table also give a 
hint about the influence of softening: as suggested above soft¬ 
ening mainly affects the spin of small subhaloes. This is not 
surprising for the reasons already indicated. 


4.2. Haloes vs. Subhaloes 

So far we have seen that the lower spin parameter distri¬ 
bution observed for subhaloes appears to be a robust result 
that does not depend upon the choice of halo finder, gravita¬ 
tional softening or mass resolution. Here we explore a possi¬ 
ble physical origin for the lower subhalo spins. 

The left and right panel of Figure IT] show the Peebles spin 
and Bullock spin distributions of subhaloes as well as haloes 
in all the simulation boxe s and the re spective fitti ng functions 
from|Bett et al.| ( |2007| i, |Bullock et al.] ( |2001| i and |Onions et al.| 
( |2013| ). To assure robustness, only structures resolved with at 
least 300 particles are included. For haloes we remove par¬ 
ticles contained within substructures when we calculate the 
spin. We make fits to the histograms in Figure [I] using Equa¬ 
tion [2] and Equation [5] The respective fitting parameters are 
given in Table [4] 

Figure [T] Table [2] Table [3] and Table [4] all show that the 
subhalo spin distribution is different from the halo spin distri¬ 
bution. This further confirms earlier results by Onions et al. 
(2013) who found that the spin distribution of subhaloes in 
the Level 4 resolution Aquarius simulation (Springel et al. 
2008) is significantly different to the one derived by Bett et al. 
( |2007[ ) for haloes in the Millennium Simulation. 

Our results reveal new information about the spin of sub¬ 
haloes. As the box size grows from BoxA to BoxD, the dis¬ 
crepancy between Ao of haloes and subhaloes decreases grad¬ 
ually, i.e. the scale of the effect is mass dependent with larger 
subhaloes tending to have higher spin. One possibility is that 
in small simulations such as BoxA or isolated halo models 
such as Aquarius-A to E studied by Onions et al. ( j2013] >. large 
substructures are generally absent. In the next section we will 
demonstrate that subhalo spins increase with subhalo mass 
while halo spins do not have a significant mass dependence. 

As an aside, it sho uld be noted that wh i le our fits do not 
mat ch those given by Bullock et al.| ( |200T| i, |Bett et al.| ( |2007| > 
and |Onions et al.| (|201 3|> exactly they are within the range of 


results covered by these works. In practice, previous work 
does not arrive at an agreement on the exact value of haloes’ 
spin. Most of these studies fix Aq in the range 0.031 — 0. 045, 
with <7 betw ee n 0.48—0.64 (see section 4.1 and Fig.7 in Shaw 
|et al.|2006) . In Bett et al. ( |2007| >, they found median values of 
A med = 0.0367 — 0.0429 for the entire population of haloes, 
depending on the definition of halo, and X me d = 0.043 for 
the catalogue of haloes they refined. Different sets of simu¬ 
lations and halo selection criteria may result in t his v ariation 
of the recovered spin parameter. As Subsection 4.1 demon¬ 
strated, such factors as mass resolution and gravitational soft- 
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Figure 1. The Peebles (upper) and Bullock (lower) spin distribution of SUB¬ 
FIND (sub-)haloes at z = 0. The colored solid histograms show the results for 
subhaloes, and black dashed histogram for haloes. Two lines show the best 
fitting functions from|Bullock et a l.| (2001) (black dashed line) and Onions 
|et al.f(201 3 1 (red solid line). Only (sub-)haloes with at least 300 particles 
have been included. For total numbers refer to Table[I| 

ening influence the spin. On the other hand, the discrepancy 
between the spin of haloes and subhaloes within our simu¬ 
lations is much larger than the bias among simulations. It 
should therefore be regarded as an intrinsic physical property 
rather than a result of different data sets. 

4.3. Mass Dependence 

As we suggest above, mass dependence can explain the dis¬ 
crepancy between the spins of haloes and subhaloes. To val¬ 
idate this, we further explore the mass dependence of (sub- 
)halo spin. Figure [2] and Figure R] show the two dimensional 
histogram of spin against (sub-)halo mass. They present a 
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Table 3 

Parameters recovered by SUBFIND for the spin distribution with different mass resolution. 


Peebles Spin 


BoxA (N p > 300) 
BoxA(7Vj, > 2400) 
BoxA_S2(Af p > 2400) 
BoxLo(Af p > 300) 


haloes subhaloes 

_Ao_ a _Ao_ a 

0.0371 ± 0.00021 2.59 ±0.047 0.0254 ± 0.00049 2.92 ±0.17 

0.0393 ± 0.00062 2.53 ±0.13 0.0281 ± 0.0011 2.78 ± 0.33 

0.0389 ± 0.00069 2.48 ±0.14 0.0288 ± 0.0010 2.51 ± 0.29 

0.0392 ± 0.00073 2.38 ±0.14 0.0293 ± 0.0008 2.12 ±0.23 


Bullock Spin 


BoxA(JV P > 300) 
BoxA(JV p > 2400) 
BoxA_S2(Ai p > 2400) 
BoxLo(lVp > 300) 


haloes 

_ * * 

0.0310 ±0.00024 0.629 ±0.007 

0.0324 ± 0.00060 0.626 ± 0.017 

0.0327 ± 0.00056 0.622 ± 0.016 

0.0337 ± 0.00049 0.616 ± 0.012 


subhaloes 

K _ * 

0.0167 ± 0.00019 0.664 ±0.009 

0.0183 ± 0.0008 0.738 ± 0.035 

0.0198 ± 0.0009 0.662 ± 0.040 

0.0207 ± 0.0007 0.582 ±0.028 


Table 4 

Parameters for the spin distribution recovered from SUBFIND haloes and subhaloes in different mass ranges. 


Peebles Spin 


haloes 


subhaloes 


Ao 


Ao 


Bett et al. < 2C 

07 


0.04326 ± 0.000020 

2.509 ±0.0033* 

- 

- 

Onions et al. 

2013 

- 

- 

0.028 

3.64 

BoxA 

BoxB 

BoxC 

BoxD 



0.0371 ± 0.00021 
0.0384 ± 0.00020 
0.0404 ± 0.00018 
0.0411 ± 0.00014 

2.59 ± 0.047 
2.63 ±0.043 
2.57 ± 0.037 
2.47 ± 0.027 

0.0254 ± 0.00049 
0.0285 ± 0.00032 
0.0308 ± 0.00033 
0.0346 ± 0.00030 

2.92 ± 0.17 
2.68 ± 0.09 
2.75 ± 0.09 
2.54 ±0.07 

Bullock Spin 



haloes 

*n 

a 

subhaloes 

A' 

Bullock et al 

(2001 

0.035 ±0.005 

0.5 ±0.3** 

- 

- 

Onions et al. 

(2013 

- 

- 

0.018*** 

0.70 

BoxA 

BoxB 



0.0310 ± 0.00024 
0.0332 ± 0.00025 

0.629 ±0.007 
0.632 ± 0.006 

0.0167 ± 0.00019 
0.0199 ± 0.00020 

0.664 ± 0.009 
0.669 ± 0.008 


BoxC 

BoxD 


0.0367 ± 0.00028 
0.0393 ± 0.00027 


0.622 ± 0.006 0.0232 ± 0.00023 0.645 ± 0.008 
0.604 ±0.006 0.0283 ± 0.00025 0.610 ± 0.007 


In Bett et ah 12007), the parameters have much smaller uncertainty because there are much lager popula¬ 
tions of halos in the work (17709121 raw FOF halos including 1332239 ’’clean” ones). 

* In Bullock et al. 1 2001 1 , the parameters have larger uncertainty because they use less haloes for fitting 

(only 5UU haloes). 

’** Note that these parameters differ from the ones originally stated in Onions et al. (2013). The original 
values were derived using an incorrect fitting routine. The values stated here are the correct values fitted to 
the Aquarius L4 data set. 


straightforward picture of how the spin distribution changes 
with mass. We use four simulations to expand the mass range. 
Contours for each simulation at the same redshift are normal¬ 
ized and stacked together into one subplot. We then divide 
the sample into 40 bins by Zog(sub-)halo mass, fit the distri¬ 
bution by Equation [2] and Equation [5] (if the sample volume 
in that bin is large enough). Ao(Aq) in each mass bin is cal¬ 
culated and marked on the plots with a cross, and a linear 
fit to Ao(Ag) against mass is indicated by the red solid line. 
From the subplots corresponding to redshift 0 (top left and 
top right), we can see that, for subhaloes, Ao(Ag) clearly in¬ 
creases with increasing subhalo mass. In contrast, the spin 
distribution of haloes is almost independent of mass. The in¬ 
creasing subhalo spin with mass is even more pronounced re¬ 
lationship for the Bullock spin parameter shown in Figure [3] 
because for the Bullock spin the (sub)halo mass has a higher 
weight (c.f. Equation TjEquation [4]). This results in a larger 
discrepancy between haloes and subhaloes at the low mass 
end. In Onions et al. < 2013] ), their samples are from a Milky 
Way like re-simulation, which contains subhaloes similar to 
those found in BoxA. So the subhalo spin distribution in their 


BoxA. [Onions et ak (2013) sug- 
anism that drives tills'dil 


work is closer to that from 
gest that the physical mechanism that drives this difference is 
mass stripping. Subhalo particles with high angular momen¬ 
tum are stripped preferentially which leads to a decrease in 
the spin parameter. Subhaloes with low mass are usually the 
ones stripped most severely. Our results strongly support the 
claims [Onions et ak[ ( j2013[ ) made. 

The slightly positive slope of A 0 , h aio{M ha i o ) and 
Aq haio(Mhaio) is inconsistent with some previous work. 
They f ound that Xo.haio is constant or has a slightly negative 
slope ((|Bett et al.|2007 Maccio et al.|2007| c.f.)). However we 
should not forget this trend includes th e effect of syst ematic 
bias between simulations. Figure 10 in Maccio et al. (2008) 
shows that X' 0 (Mh a io) has a slope of 0.005 in a simulation us¬ 
ing WMAP1, which is the same cosmology as that used here. 


4.4. Radial Dependence 

To understand further whether tidal stripping of high angu¬ 
lar momentum material could cause the lower subhalo spin 
distribution, we investigate the radial dependence of the sub¬ 
halo spins. Subhaloes located closer to the center of their 
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Figure 2. 2 dimensional histogram of Peebles spin against SUBFIND (sub-)haloes mass. Contours with different colors represent different level of number density 
(calculated as dN/dlog\/dlog(M/10 10 h~ 1 Mq)/10000) as indicated in the color bar above. Subplots in left column are statistic for haloes and right column 
for subhaloes. Two plots in a same row are from the same snapshot. The redshift of each row increase from top to bottom respectively. There are in fact four 
parts of contour in each subplot, which comes from our four simulations respectively. The cross scatters represent Ao in every mass bin. The red thick lines are 
linear fitting to Ao against mass. The black dashed line indicates a position of log A = — 1.5(A ~ 0.032) as a standard for comparison. 
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Figure 3. 2 dimensional histogram of Bullock spin against SUBFIND (sub-)haloes mass. Contours with different colors represent different level of number 
density (calculated as dN/dlogX'/dlog(M/10 10 h~ 1 Mq)/10000) as indicated in the color bar above. Subplots in left column are statistic for haloes and right 
column for subhaloes. Two plots in a same row are from the same snapshot. The redshift of each row increase from top to bottom respectively. There are in fact 
four parts of contour in each subplot, which comes from our four simulations respectively. The cross scatters represent Aq in every mass bin. The red thick lines 
are linear fitting to Aq against mass. The black dashed line indicates a position of logX' = — 1.5(A ; « 0.032) as a standard for comparison. 
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host halo are likely to have undergone stronger tidal stripping 
than those nearer the virial radius. |Onions et al.| ( |2013) has 
done some tests to support their argument, e.g. they analyze 
the average spin parameter of subhaloes at different distances 
from the center of the host halo. Here we perform a more 
detailed test. We stack subhalo samples from four simula¬ 
tions together and then make the two dimensional histogram 
of spin against their centric distance. Then in each radial 
bin we fit the subhalo sample using Equation [2] and Equa- 
tionBl Finally we make a linear fit as Ao(r) = cr/R V i r + d or 
Ao(r) = cr/R V i r + d. The results are displayed in Figure [5] 
and show that the spin of subhaloes is suppressed close to the 
center of the host halo. This is consistent with the argument 
that subhaloes loose their high angular momentum particles as 
they are stripped of their outer layers after infall into a main 
halo. 


4.5. Redshift Dependence 

So far all our analysis was conducted on the 2 = 0 snapshot. 
However, the spin distribution is known to change with red- 
shift (Hetznecker & Burkert| |2006[ ). Hence, we investigated 
the redshift dependence of the halo vs. subhalo spin distribu¬ 
tion offset. In Figure [2] and Figure p] we give the Peebles and 
Bullock spin distributions of (sub-)haloes at z = 0, z = 0.84, 
z = 2.38 and z = 5.21, respectively. The spin of subhaloes 
at the lower mass end decreases significantly with time, while 
the spin of massive haloes increases slightly with time. We 
calculate the Ao(Ag) in each mass bin and then use a linear 
function log\o = a*log(M/10 10 h~ 1 Mo)+b(\ l 0 for Bullock 
spin) to fit Ao(Ag) against mass. We list the parameters of each 
fitting line in Table [5] It is clear the difference between the 
halo and subhalo spin distribution increases with time. This 
is consistent with the argument that tidal stripping causes the 
difference. Affected by stripping, subhaloes loose more and 
more high angular momentum particles as time passes. 

To confirm this result we checked how the spin of a sin¬ 
gle halo changes with time. We randomly select 6 subhaloes 
and plot their spin against redshift. We constrain the samples 
selected so that they are likely to be the subhaloes heavily 
stripped. As expected, heavily stripped subhaloes should have 
low mass (here we choose masses less than 10 9 H~ 1 Mq) and 
have long histories, forming prior to redshift 8. As shown in 
Figure [5] the spin of subhaloes declines at low redshift. We 
have checked many more subhaloes not displayed in Figure [5] 
and we find that most of them display the same trend. We 
also plot a dashed line for a halo as reference. The halo is in 
the same mass range as those subhaloes selected. Its spin al¬ 
most does not change at low redshift. This piece of evidence 
strongly supports the claim that subhaloes suffer from strip¬ 
ping, loosing their spin over time. 

The information at the high redshift end in Figure [5] is not 
reliable since the progenitors don’t contain very many parti¬ 
cles. This results in the large fluctuations seen here. 


5. DISCUSSION & CONCLUSIONS 

In this work we compared the spin distribution function of 
haloes and subhaloes in sets of cosmological box simulations. 
We found that the halo spin distribution function is well fit¬ 
ted by the parametrisations given by Bullock et al. (2001) and 
|Bett et aL] ( |2007) for the Bullock and Peebles ( |Pee ~les|1969[ i 
spin parameter, respectively. For the subhalo spin distribu¬ 
tion function, however, the typical spin of a small subhalo is 
significantly lower. This was previously suggested by|Onionsj 



1.8 


1.5 


1.2 


0.9 


0.6 


0.3 




- log\'„ =0.17r/U„ ir -1.80 


-0.5 


O 



1.8 


1.5 


1.2 


0.9 


0.6 


0.3 


Figure 4. The average Peebles (upper) and Bullock (lower) spin distribution 
of SUBFIND subhaloes against their relative distance from the host halo cen¬ 
ter. r is the distance of subhaloes from the center of their host halo. R v i r 
is the virial radius of host halo. The crosses indicate the peak of every ra¬ 
dial bin and the red line is a linear fitting to these marks. Contours with 
different colors represent different level of number density ( calculated as 
dN/dlogX(orX')/dlog(r/R v i r )/10000 ) as indicated in the color bars to 
the right. 


et al. (]2013j) for the subhaloes within Milky Way-like haloes, 
but is confirmed here for a full cosmological volume. 

We investigated the origin of the difference between the 
halo and subhalo spin distributions. We examined the influ¬ 
ence of (sub-)halo finder, spin parametrisation and resolution 
to confirm that these factors are not the origin of the difference 
between the halo and subhalo spin distributions. In this pro¬ 
cess we confirmed the difficulties that the AHF finder has in 
recovering substructures reliably where the density contrast 
between the main halo and the subhalo is low. We recom- 
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Figure 5. The evolution of Peebles (upper) and Bullock (lower) spin of six 
selected subhaloes from BoxA. Subhaloes have at least 300 particles and are 
less massive than \t)‘ J h l VL.,. All selected subhaloes form prior to redshift 
8. Solid colored lines represent different subhaloes. The black dashed line 
shows the spin of a halo changing with time as a reference. 

mend that the AHF halo finder should be treated with caution 
in situations where a complete unbiased sample of the subhalo 
population is required. 

In this paper, we have argued that the difference between 
the spin distributions is physical and it is caused by tidal strip¬ 
ping of subhaloes removing high angular momentum mate¬ 
rial. This argument is strongly supported by three pieces of 
evidence presented here. Firstly, subhaloes tend to have lower 
spin when compared to haloes of the same mass. This dis¬ 
crepancy gets larger towards the low mass end. Secondly, the 
spin distribution of subhaloes is radially dependent within a 
host halo. Subhaloes closer to the host halo center, which are 


Table 5 

Parameters for fitted lines for spin against (sub-)halo mass 




log Ao 

= a*log(M/10 10 h 

~ 1 M Q ) + b 




Halo 


Subhalo 



a 

b 

a 

b 


2 = 0 

0.013 

-1.41 

0.039 

-1.53 

Peebles spin 

2 = 0.84 

0.0015 

-1.42 

0.037 

-1.49 

2 = 2.38 

-0.010 

-1.43 

0.018 

-1.44 


z = 5.21 

-0.017 

-1.46 

0.018 

-1.41 



log\' n 

= a*log(M/10 10 h 

~ 1 M Q ) + b 




Halo 


Subhalo 



a 

b 

a 

b 


2 = 0 

0.028 

-1.47 

0.069 

-1.68 

Bullock spin 

2 = 0.84 

0.021 

-1.46 

0.066 

-1.61 

2 = 2.38 

0.0073 

-1.42 

0.047 

-1.47 


z = 5.21 

-0.0068 

-1.39 

0.037 

-1.37 


expected to have been more tidally stripped, have lower spin 
than those closer to the virial radius. Thirdly the difference 
between halo and subhalo spin increases with time and hence 
is being caused by a dynamical effect, such as tidal stripping. 

In summary, we have demonstrated that subhaloes typically 
have lower spin than haloes because tidal stripping removes 
their highest angular momentum material. This can have an 
important consequence for galaxy properties that require spin 
parameter information. Galaxy properties that are related to 
spin are more likely to be correlated to the spin of the subhalo 
before infall and not necessarily to its present value. 
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